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ABSTRACT 


The Phase-Induced Amplitude Apodization Complex Mask Coronagraph (PIAACMC) provides an efficient way 
to control diffraction propagation effects caused by the central obstruction/segmented mirrors of the telescope. 
PIAACMC can be optimized in a way that takes into account both chromatic diffraction effects caused by the 
telescope obstructed aperture and tip/tilt sensitivity of the coronagraph. As a result, unlike classic PIAA, the 
PIAACMC mirror shapes are often slightly asymmetric even for an on-axis configuration and require more care 
in calculating off-axis shapes when an off-axis configuration is preferred. 


A method to design off-axis PIAA mirror shapes given an on-axis mirror design is presented. The algorithm 
is based on geometrical ray tracing and is able to calculate off-axis PIAA mirror shapes for an arbitrary geometry 
of the input and output beams. The method is demonstrated using the third generation PLIAACMC design for 
WFIRST-AFTA telescope. Geometrical optics design issues related to the off-axis diffraction propagation effects 
are also discussed. 
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1. INTRODUCTION 


The WFIRST-AFTA coronagraph is designed to have powerful high contrast performance for direct imaging of 
exoplanets by their reflected light. Any coronagraph system considered for the WFIRST-AFTA project should 
solve the main design problem how to remove the star light diffracted from the telescope central obstruction. The 
Phase-Induced Amplitude Apodization Complex Mask Coronagraph (PIAACMC) coronagraph! based on lossless 
pupil apodization by the beam shaping concept? (PIAA coronagraph) combined with a complex amplitude focal 
plane mask gives a solution of this problem. 


The PIAACMC coronagraph can be optimized in a way that takes into account both chromatic diffraction 
effects caused by the telescope obstructed aperture, and tip/tilt sensitivity of the coronagraph,? and still attain 
the performance required by AFTA.* The typical design process for PIAA mirrors is to first optimize their 
shapes for an abstract (unphysical) on-axis collimated space configuration and then convert this design to a real 
unobstructed off-axis physical configuration. The conversion to off-axis is not trivial and the goal of this paper is 
to present a method to calculate PIAA optics for an arbitrary off-axis geometry of the optical system based on the 
optimized on-axis optics shapes. By ” arbitrary” , we mean different off-axis geometries for the system entrance and 
exit beams as well as a possible beam compression that the system should perform. Though diffraction critically 
affects the final system performance, our consideration will be based on geometrical optics only, assuming that the 
diffraction propagation effects are taken into account during the on-axis system optimization and are unaffected 
by conversion to off-axis. 


collimated output beam 












F (f sin B, d-f cos B) PIAA M2 


(0.0) 








esgrrcccteryrn eee sree (0, d) 





PIAA M1 








collimated input beam FLAT 


F — the main PIAA focus. 

d — the distance between PIAA mirrors. 
a — the incoming beam angle. 

B — the outgoing beam angle. 


The collimated input/output beams are used for remapping calculations. 


Figure 1. Optical layout. 


2. OPTICAL LAYOUT 


The proposed optical layout? of the PIAA system for the WFIRST-AFTA PIAACMC coronagraph (Fig. 1) 
includes two PIAA mirrors M1 and M2. In this collimated beam-to focus system the collimated input beam is 
reshaped by the M1 mirror and is focused by the M2 mirror. The angle of the input beam a is equal to 20° 
(approx.) and the output beam angle ( is equal to 25° (approx.). The input beam diameter is about 40 mm, 
the distance d between mirrors is about 750 mm, and the focal length f of M2 (the distance between the M2 
center and the system focus) is about 1200 mm. The M1 mirror also performes the compression of the input 
beam needed to match the focal plane PSF size with the diameter of the focal plane occulter. As a result the 
beam diameter on the M2 is approximately equal to 15 mm. 


We define a coordinate reference frame shown in Fig. 1 with the origin at the center of PIAA M1 and the z 
axis parallel to the central ray between PIAA M1 and M2. All mirror shapes will be expressed in this reference 
frame as functions of (a, y). 


The base shape of the M1 and M2 mirrors (i.e. without any PIAA remapping) are a plane and an off-axis 
parabola (OAP) respectively. All the following mirror shapes are calculated as corrections to these shapes that 
are given by: 


PLANE(a,y) = y x tan(a/2), 
f? sin? 6 —2? —(y — f sing)? 


OAKEY) S 2F (1 + cos B) 


+d. (1) 
The base shapes are computed such that they focus the collimated input beam to the desired system focus whose 
coordinates are (x,y,z) = (0, fsinG,d — f cos) as shown in Fig. 1. Such a choice of base shapes turns the 
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Figure 2. Extended remapping function R2(r1) and inverse remapping function Ri(r2): the working beam radii rp; and 
rpg are shown. 


calculation of these shapes into a trivial procedure that is independent on a particular focus position and the 
beam compression/extending rate. 


Two additional mirrors are used to perform remapping calculations. They are the flat mirror FLAT and the 
off-axis parabola OAPy: 
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—f? sin? B+ 2? + (y — f sin 8)? 
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FLAT(az,y) = 
OAP;(x,y) = 


Both mirrors form collimated input and output beams that are parallel to the optical axis between PIAA M1 
and M2 mirrors. 


3. OFF-AXIS MIRROR SHAPES 


In this section the procedure to calculate off axis PIAA shapes is given. The shape of both mirrors can be 
presented as a sum of three terms: 


Mi (a2,y) = PLANE(a,y) + fi(r) + OAT; (a, y) 
Mo(a,y) = OAP(a2, y) + fo(r) + OAT2 (a, y). (3) 


In Eq. 3 PLANE(a, y) and OAP(a, y) are base shapes of PIAA M1 and M2 mirrors (Eq. 1), fi(r) and fo(r) 
are the circularly symmetric mirror terms that are mainly responsible for the beam remapping in the on-axis 
collimated beam-to-focus system, and OAT; (x,y) and OAT»(x, y) are the off-axis PIAA correction terms that 
preserve the correct PIAA apodization when mirrors are converted to off-axis. 


3.1 Remapping functions 


The ray tracing calculations of the PIAA optics are based on the remapping function concept. For an on-axis 
source the remapping function establishes a correspondence between a ray position ry = (21, y1) in the entrance 
pupil and the position rg = (x2, y2) of the same ray in the exit pupil: 


x2 = %2(%1, 41) 


yo = y2(1,y1)- (4) 


Note, that the origin of the reference frames in Eq. 4 matches with the position of the central ray in the entrance 
and exit pupils. 

Since we take the on-axis mirror shapes as given, the ray positions in the entrance and exit pupils, and thus 
2D remapping functions x2(21, y1) and y2(x1, y1), can be readily calculated with a ray tracing procedure. 


The starting point to design mirror shapes for the WFIRST-AFTA coronagraph was the diffraction-based 
optimized solution obtained for an on-axis collimated beam to collimated beam PIAA system without beam 
compression.® Although this solution is not circularly symmetric because of the central obstruction asymmetry, 
the deviation from the symmetry for both mirrors is small. That makes it possible to define the azimuthally 
averaged remapping function 


ro(r1) =< 4/23(r1 cos ¢, 1 sin d) + y3(ri cos ¢, 71 sin ¢) >¢, (5) 


that determines the radial mirror terms f;(r) and fo(r). In Eq. 5 <... >4 means azimuthal averaging. 


The remapping function r2(r1) constructed in according with Eq. 5 provides a good initial approximation 
to calculate the off-axis mirror shapes for all the rays within the working beam. However, the discontinuity of 
derivatives of the remapping function at the beam edge and an absence of the remapped values outside of the 
working beam create numerous difficulties for the ray tracing calculations that suppose to work with real square 
arrays of points. 

To avoid these problems we created an extended remapping function Ra(r1) (Fig. 2) by transposing the 
remapping function r2(r1) onto the intervals [rb;, R] and [rbz, R] backwards, beyond the working beam radii rp1 
(on the PIAA M1) and rp (on the PIAA M2), and making sure that the derivatives match on the beam edge. 
The Ro(r1) is given as follows: 





Ro(r1) = r2(71) for T. <To1; (6) 
ry =ro t+ (lh tairi)ry 
ro = 72 + (1+ aars)r4 for ro. <1 < R, 





rpg — Ty = Ro(ro1 — 14) 


where 
Bite (R= 101)/To1 — Es 
Tol 
= aA 
ee (R — ro2)/Toe 
Tb2 


and R is the radius of an area big enough to cover all arrays involved in the mirror shapes calculation. 

It should be noted that for a system with the beam compression rp2 is not equal to rp;. The desired beam 
compression can be also included in the final expressions for 2D remapping functions X2(x1, yi) and Y2(#1, y1) 
through the rescaling factor rp2/rp1: 


rT 
Xo(%1,y1) = Fate, 91) 
T 
Yo(z1,y1) = = ua(tiy); (7) 


where x2(21,y1) and yo(x1, yi) are remapping functions for the system without the beam compression. 
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Figure 3. Mirror shapes for WFIRST-AFTA coronagraph: radial mirror terms, the focal distance f is equal to the distance 
between mirrors d. 


3.2 Optics shapes calculation 


The following iteration procedure is proposed to calculate off-axis PIAA mirror shapes given on-axis ones: 


1. Calculate the 2-D remapping functions x2(21, y1) and y2(x1, yi) by ray tracing on-axis mirror shapes and 
azimuthally average them to get the radial remapping function R(r1). 


2. Calculate the 2-D remapping functions X2(x1, yi) and Y2(x1,y1) for the compressed beam. 
3. Determine base shapes PLAN E(x,y) and OAP(a, y) of mirrors per Eq. 1. 


4. Calculate radial mirror profiles f(r) and f2(r) for on-axis collimated beam-to-focus system from R2(r1) 
(Section 3.3). 


5. Start a loop. For i-th iteration step construct an estimate OAT, ;(x,y) as sum of small correcting terms 
dM, ,;(x, y) calculated during previous iteration steps (Section 3.4): 


OAT;,i(2,y) = 5 dM15(2,y), 
j=0 
OAT, o(x,y) = 0. (8) 


6. Determine off-axis PIAA term OAT>,;(x, y) that provides equal optical path length (OPL) for all the rays 
propagated from the source to the main PIAA focus system. This step is performed assuming the off-axis 
PIAA term OAT;(z, y) is equal to OAT}; (x, y). 


7. Close the loop by returning to step 5. 


3.3 Radial mirror terms 


For the case of a radially symmetric collimated-beam to collimated-beam PIAA system, the mirror shapes can 
be directly computed from R2(r1) by numerically integrating a pair of differential equations.” 


In this section we consider the case of the on-axis collimated-beam-to-focus radially symmetric system. The 
proposed solution is similar to the procedure we use to calculate 2D off-axis later in this paper. 


The radial mirror profiles f,(r) can be derived from 


" dfi(t 
iy =10)+ [| Sa, (9) 
0 t 

if the derivative dfi(r)/dr is known. The derivative dfi(r)/dr can be easily estimated for a collimated beam-to- 
focus system by using the ”desired” normal (i.e. one that generates the desired remapping in according with 
Eq. 5) to the M1 surface. 


In according with the reflection law the ”desired” normal nj = (n’.,n‘,) to the surface of the M1 mirror can 


be written as ‘ - 
1s 12 
= 





ni(t1) = (n,,n,) = (10) 
Tl1s T12 

where ry; is the vector from a point ry on the M1 surface to the source, and the vector ry2g connects point rz 
with the remapped point rg. The position of the remapped point rg on the M2 mirror is determined by the 


remapping function R2(r1). The derivative df(r)/dr is given by 


df1(r) = _ Pr (11) 


a 7 
dr nt 





To derive both fi(r) and f2(r) radial profiles an iterative approach is used. In this approach, the ” desired” 
normal n; is estimated from the previous approximation for the M1 and M2 shapes and the M1 shape is obtained 
by solving Eq. 9. The next approximation of the M2 shape is calculated assuming that for all the rays propagated 
from the source to the main PIAA focus OPL should be equal to the OPL of the principal ray (Fermat’s principle). 
The OPL of the principal ray is known from the initial conditions. Radial mirror terms for the WFIRST-AFTA 
coronagraph are shown in Fig. 3. 


3.4 Calculation of dM,(z, y) 


To calculate the corrective term dM1,;(x, y) (Section 3.2) we should: 


1. For given 2D remapping functions (Eq. 7) determine the ”desired” normal nj = (ni,,nj,,n4) to the M1 
surface assuming that M1 and M2 mirror shapes match the current shapes approximation. Similarly 
Eq. 10, the ” desired” normal 


ni(x, y) a (tne) = 





eee (12) 
Tls T12 
where vectors ri, and ry2 are determined in Section 3.3. To determine a particular ray position on the 
M2 surface the geometrical propagation in direction from the exit pupil to the off-axis parabola OAPy 
and then to the M2 mirror should be used. The ray position in the exit pupil can be determined as the 
remapped (in accordance with Eq. 7) ray position in the entrance pupil. By the entrance pupil we mean 
a plane upstream of the flat mirror FLAT(x,y) (Eq. 2). The exit pupil is the plane downstream of the 
off-axis parabola OAP;(z, y) (Eq. 2) optically conjugated with the M2 mirror. Both the entrance and exit 
pupils are orthogonal to the collimated input/output beams and all the rays are considered to be collimated 
between the entrance pupil and the flat mirror and between the exit pupil and OAPy. 
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2. Obtain the ”current” normal (i.e. gradient vector of the surface) ne = (no, nf, 
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Figure 4. PIAA M1 (top) and M2 (bottom) off-axis shapes for the WFIRST-AFTA coronagraph. 
Working beam radii ry: and ry2 are shown on crosscuts on the right. 


3. Calculate derivatives 0dM1,;(x, y)/Ox and 0dM, ; (x, y)/Oy: 


ddMy (x, y)/Ox = —ni,/ni, + n§/ng, (13) 
ddMy,s(x, y)/Oy = —ni,/ni, + n§/n§. (14) 


4. Calculate the dM1,;(x,y) through a line integral 


(x,y) 
datsesy) = [OED ga 5 ay (15) 


Since, for any physically possible remapping the integration in Eq. 15 along any path that connects the 
mirror center with the point (x,y) should give the same result, dM1,;(x, y) values can be also averaged for 
different integration paths. 


5. To reduce the M1 errors caused by the numerical integration/interpolation an appropriate smoothing 
algorithm should be applied for each dM, ;(x, y) estimation. 


3.5 Calculation of OAT»(z, y) 


To determine the OAT) ;(x, y) we use an iterative procedure II that corrects the OPL of the output rays on every 
step of the main iteration. Note that iteration procedure II is different from the main iteration discussed earlier 


in Sec. 3.2. 
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Figure 5. The output beam geometrical OPD(x,y) for the on-axis optimized shapes in projection on the M1 (top) and the 
PHASE(x,y) mirror term (bottom). Note that OPD(ax,y) and PHASE(z,y) have different scales. 


During k-th step of the iteration II the function dM2,4(x2, y2) is calculated: 


OP Eas, in) SOPLO,0 
dMo,4(2, 42) = oe a (16) 


where OPL;(z, y) is the optical path length of the ray that starts at the source position, reflects from the current 
approximation of the M1 shape at the point (x1, y1), reflects from the current approximation M2 shape at the 
point (v2, y2) and then hits the system focus. Similary to OAT) ,;(z,y), the off-axis term OAT2,(x, y) can be 
determined as 


OAT? (£2, y2) = S dM2,4(x2, 42), (17) 
k 


OAT? 0(x2, y2) = 0. 


The iteration stops as soon as the OPD for all the rays that cross the M2 surface becomes negligible. Because 
of the remapping, the sampling of OAT2(x2, y2) on the M2 surface appears to be inhomogenious, even if the 
original arrays sampling on the M1 surface is homogenious. To get homogeneous OAT>(x2, y2) sampling on the 
M2 surface a triangulation procedure should be used. 


4. MIRROR SHAPES 


The calculated mirror shapes for the WFIRST-AFTA coronagraph are shown in Fig. 4. To present the shapes 
deviation from common optical shapes, the parabolic terms given by 


P,(x, y) = 0.171215209 x y + 0.24227703 x (x? + y”),  subtructed from PIAA M1 


Po(x, y) = 0.228083501 x y + 0.64898643 x (2 + y2)  subtructed from PIAA M2. (18) 


were pre-subtracted from the mirror surfaces. Note that these parabolas are not the base optics shapes considered 
earlier and are not chosen based on their optical properties, they were chosen simply to highlight higher order 
departures visually. 


An additional point to emphasize is that the diffraction propagation introduces radial oscillations both in 
amplitude and in phase of the output beam (Fig. 5,6) which the initially optimized on-axis mirror shapes take 
into account and correct for. 


As a result the geometrical OPD of the output beam is not equal to zero according to a purely geometrical 
optics model (Fig. 5), although it appears to be equal to zero if we account for diffraction propagation. To take 
into account the above effect an additional term for the M2 mirror shape PH AS E(x, y) should be introduced: 


Ma(a,y) = OAP(a,y) + fo(r) + OATo(a, y) + PH ASE(z, y), (19) 


where OAT»(x,y) term describes the M2 shape assuming that the wavefront after M2 mirror is spherical, and 
PHASE(x,y) takes into account the OPD change caused by the diffraction propagation in the system. To calculate 
the M2 shape in this case the equation Eq. 16 should be modified in the following way 


_ OPL; (x2, y2) — OPL(0,0) — OPDa(z, y) 


dMo2,x (2, y2) = 5 ‘ (20) 


where OPDa(z, y) is the geometrical OPD in the output beam for the baseline on-axis system. 


The second problem related to the diffraction caused oscillations in the mirror shapes is the dramatic slow 
down of the procedure to calculate M2 shape when the residual OPD of the output beam reaches 107!° m level. 


5. RAY TRACING CHECK 


The derived optical shapes with the ”phase” term PHASE(x,y) = 0 have been checked with the ray tracing 
code that has been developed to simulate the previous generations of PIAA optics manufactured by AXSYS 
Inc.® and Tinsley Inc.® Results of the ray tracing simulation shown in Fig. 6 


The complex amplitude in the exit pupil of the PIAA has been obtained by analysing the distribution of 
the 2 x 108 random rays is enough to estimate the amplitude profile with rms of residuals of about 0.09. The 
derived amplitude profile clearly shows the radial amplitude oscillations. These oscillations would be corrected 
if the diffraction propagation is considered. 


The phase deviations in the exit pupil are only limited by the residual OPD of the output beam caused by 
the convergence of the M2 calculation algorithm discussed above. Their rms reaches the /3000 level that is 
comparable with the physically achievable limit for the optical quality. 


6. SUMMARY 


In this paper we presented a method to calculate off-axis mirror shapes for PIAA coronagraph assuming that 
the on-axis mirror shapes are known. The method is based on geometrical ray tracing and is able to calculate 
off-axis PIAA mirror shapes for an arbitrary geometry of the input and output beams. 
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Figure 6. Three dimensional ray tracing check of the PIAA mirror shapes with the term PH AS E(a, y)=0: 
a) the exit pupil amplitude and phase; 
b) the amplitude and phase crosscuts through the exit pupil center. 
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